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The jamming and percolation for two generalized models of random sequential adsorption (RSA) 
of linear fc-mers (particles occupying k adjacent sites) on a square lattice are studied by means 
of Monte Carlo simulation. The classical random sequential adsorption (RSA) model assumes the 
absence of overlapping of the new incoming particle with the previously deposited ones. The first 
model LKd is a generalized variant of the RSA model for both fc-mers and a lattice with defects. 

Some of the occupying fc adjacent sites are considered as insulating and some of the lattice sites are 
occupied by defects (impurities). For this model even a small concentration of defects can inhibit 
percolation for relatively long fc-mers. The second model is the cooperative sequential adsorption 
(CSA) one, where, for each new fc-mer, only a restricted number of lateral contacts 2 with previously 
deposited fc-mers is allowed. Deposition occurs in the case when z < (1 — d)zm where Zm = 2(fc -|- 1) 
is the maximum numbers of the contacts of fc-mer, and d is the fraction of forbidden NN contacts. 

Percolation is observed only at some interval kmin < k < kmax where the values kmin and kmax 
depend upon the fraction of forbidden contacts d. The value kmax decreases as d increases. A 
logarithmic dependence of the type \og{km,ax) = a + bd, where a = —4.03 ± 0.22, b = 4.93 ± 0.57, is 
obtained. 
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I. INTRODUCTION 

The model of random sequential adsorption (RSA) is 
very popular in studies of the colloidal deposition of par¬ 
ticles (proteins, nanoparticles, viruse^ biological cells 
etc.) on different types of substrate [l|. The particles 
are randomly deposited on the substrate with the pro¬ 
cess being fully irreversible, without subsequent detach¬ 
ment or diffusion. The classical variant of this model 
assumes an absence of overlap of the new incoming par¬ 
ticle with the previously deposited ones. With a large 
enough concentration of the deposited objects, they can 
form a spanning path between the opposite sides of the 
substrate and this concentration corresponds to the per¬ 
colation threshold 0. Finally, the jamming limit will be 
reached beyond which no more objects can be adsorbed. 

Very often in the simulations, the substrate is consid¬ 
ered as a discrete space, e.g. a square lattice. Great 
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effort has been devoted to studies regarding percolation 
and jamming for the RSA deposition of elongated par- 
tides, e.g., sticks Q, line segments 0, rigid rods, nee¬ 
dles [^or linear fc-mers (particles occupying fc adjacent 
sites) [1,0 ■ For clarity, hereinafter we shall use the term 
fc-mer. Problems with completely disordered fc-mers 
i&113 or partially aligned fc-mers @,[3, [HI 113 have been 
reported over the last two decades. Intensive studies have 
shown that the jamming concentration continuously de¬ 
creases as the length of the fc-mer increases On the 
other hand, the percolation threshold initially decreases 
and then increases with increasing value of fc Q- For 
a completely disordered system, a conjecture has been 
offered that percolation is impossible when fc exceeds ap¬ 
proximately 1.2x10^ Q. Direct verification of the conjec¬ 
ture is very time-consuming and problematic even with 
a high-performance computer. 

Different variants of the more general RSA models, 
taking account of the heterogeneity of substrates, inter¬ 
actions between the deposited particles and the possi- 
bility of surface diffusion have been proposed |13l - ll7j| . 
These models are more realistic in their description of 
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the experimental results for colloid particle adsorption 
on substrates characterized by a wide spectrum of bind¬ 
ing energies. 

Very often, the real surfaces are chemically heteroge¬ 
neous and contain defects [l^, moreover, the substrates 
may be prepatterned [l^. The structure of the elon¬ 
gated particles, e.g., carbon nanotubes, adsorbed on the 
substrate may also be highly heterogeneous, e.g. due to 
their chemical functionalization [^ . The jamming and 
percolation of fc-mers on disordered (or heterogeneous) 
substrates with defects, or fc-mers with defects, have also 
attracted great attention [l^ . [2ll - l^ . Two models of non¬ 
ideal lattices and objects were analyzed [s^. In the first 
model, it is assumed that the initial square lattice is non¬ 
ideal and that some fraction of the sites, d, is occupied 
by point defects (impurities). The lattice sides occupied 
by these point defects are forbidden for deposition of the 
objects. In the second model, it was assumed that the 
square lattice is perfect and some fraction of the sites in 
the fc-mers, d, consists of defects, i.e., is non-conducting. 
For both models, above some critical concentrations of 
defects, dm, the percolation is blocked even at the jam¬ 
ming concentration of A:-mers. The estimations predicted 
the absence of percolation even for non-defective systems 
when fc ^ 6 X 10^ [s^. Integration of the above models 
into one generalized model in which the defects are pre¬ 
sented both on the lattice and inside the /c-mers looks 
very attractive and promising, because in the real world, 
both the deposited objects and the substrates are, as a 
rule, non-ideal. We shall denote this model as the 
model. 

Another generalized RSA model of cooperative sequen¬ 
tial adsorption is that in which the adsorption probabili¬ 
ties are dependent on the local environment [l| is of spe¬ 
cial interest. This model takes account of the presence 
of very strong near-neighbor (NN) lateral repulsive in¬ 
teractions. In the simplest case, the constraint assumes 
that all NN locations are empty. For the deposition of 
monomers on a square lattice with complete NN exclu¬ 
sion, jamming is observed at a coverage pj = 0.3641. 
However, for this problem, percolation never actually oc¬ 
curs due to jamming. In the general case, a restricted 
number of lateral contacts with previously deposited par¬ 
ticles is allowed. In this model with partial NN exclusion 
(the Cd model), the fraction of forbidden NN contacts, d, 
may be identified with the fraction of defects that influ¬ 
ence the process of deposition. Note that the defect-free 
variant of this model {d = 0) corresponds to the classical 
RSA model. 

The goal of the present research is to investigate the 
jamming and percolation for the LK^ and Cd models of 
the deposition of linear fc-mers on a square lattice. In the 
LKd model, the simultaneous presence of defects both on 
the substrate and in the deposited objects is assumed. 
For this model, the defects on the substrate influence 
the deposition process and the defects on the fc-mers af¬ 
fect the connectivity of the system. In the Cd model, 
restricted number of lateral contacts with previously de¬ 


posited /c-mers is allowed for each new fc-mer. For this 
model, the fraction of forbidden NN contacts (or concen¬ 
tration of defects d) influences the deposition process. 

The rest of the paper is constructed as follows. In 
Section m we describe the technical details of our simu¬ 
lations. Section m presents our main findings. In Sec¬ 
tion lYl we summarize the results and conclude the pa¬ 
per. 


II. DETAILS OF SIMULATION 

The percolation and jamming behavior of elongated 
objects on a substrate was investigated using computer 
simulation. In both the LK^ and Cd models, we consid¬ 
ered a discrete two-dimensional substrate (square lattice 
L X L sides) with periodic (toroidal) boundary condi¬ 
tions. The deposited objects were linear A:-mers (par¬ 
ticles occupying k adjacent sites). Isotropic deposition 
was simulated, i.e., the A:-mers being deposited in two 
allowed perpendicular directions (vertical or horizontal) 
with equal probabilities. To extrapolate the results of 
the simulation to the thermodynamic limit (L —>■ oo), we 
performed a scaling analysis l^. 

A. LKd model 

In the LKd model, the presence of defects with a con¬ 
centration of dk on the deposited fc-mers was assumed. 
The length of the /c-mers, k, varied from 2 to 64. The 
lattice was also imperfect, i.e. some fraction of the lattice 
sites, di was occupied by defects (impurities) (Figure [T]). 
The defects hinder the adsorption of the elongated ob¬ 
jects. 

For each given concentration of the defects on the lat¬ 
tice, di, 

• we filled the lattice with randomly distributed point 
defects, 

• then, we deposited fc-mers up to the jamming state, 
using RSA, 

• subsequently, we randomly placed defects onto the 
fc-mers; the defect concentration, dk, being defined 
as the fraction of insulating sites on the fc-mers, i.e. 
the number of such defects is dkpL"^, 

• and finally, we checked whether percolation oc¬ 
curred. 

We filled the lattice at least 1000 times and found the 
probability, R{dk), that percolation would occur at given 
concentrations of defects on the fc-mers. 

In contrast to [3l|, we chose to treat spiral clusters as 
wrapping (percolating). We checked the percolation in 
two perpendicular direction and used two criteria: there 
is percolation in both directions (criterion AND), or there 



(a) Jamming state (b) Spanning cluster 

FIG. 1. Example of a jamming state of isotropically deposited fc-mers with defects on a square lattice with defects. The lattice 
size is 1024 x 1024, fc = 8, a fragment of the lattices with 128 x 128 sites is shown. The concentration of defects on the lattice is 
0.02. The concentration of defects on the fc-mers is 0.2. Online: Horizontal fc-mers are shown in red, vertical fc-mers are shown 
in blue, fe-mers belonging to the percolation cluster are shown in the same colors but with different brightness, empty sites are 
shown in white, defects on the lattice are shown in black, defects on the fc-mers are shown in the tint of the color of the fc-mer. 
Print: gray-scale. 


is percolation at least along one direction (OR). The ab¬ 
scissa of the inflection point of the curve was treated as 
the estimation of the critical concentration of defects for 
the given lattice size. We utilized a scaling relation to 
obtain the critical concentration at the thermodynamic 
limit 

4(L)-4(oo)) 

where v^. = 4/3 is the critical exponent 0. For fc < 32, 
we used the lattice sizes L = lOOfc, 200fc, and 400fc. For 
k = 64, we used only two lattice sizes L — 100/c and 200fc 
but the scaling analyzes were performed using the two 
criteria (AND, OR). 

B. Cd model 

In the Cd model, deposition occurs for the case when 
z < {1 — d)zm where Zm = 2(/c -|- 1) is the maximum 
numbers of contacts of the fc-mer, and d is the fraction 
of forbidden NN contacts. We filled up the lattice to a 
given concentration of fc-mers, p, 1000 times and found 
the probability, R{p) , that of percolation occurring. The 
abscissa of the inflection point of the curve was treated as 
the estimation of the percolation threshold for the given 
lattice size. The filling fraction of the lattice by fc-mers 
is defined as p = iVfc/L^, where N is the number of fc- 
mers. The value p changes within [0;pj], where pj is the 
jamming concentration. 

Figure [2] demonstrates typical examples of the jam¬ 
ming states for the Cd model with restricted numbers of 
lateral contacts for fc = 8 {zm = 18) and different values 
of d. Percolation is absent at a large number of forbid¬ 


den contacts, z > 14 (d = 14/18) (d^) but is present for 
smaller number of contacts, z > 11 (d = 11/18) (Fig.[2]D). 
The decrease of lateral repulsion between the fc-mers re¬ 
sults in increasing connectivity of the system and in the 
appearance of a percolation cluster (Fig. [2j)). 

We used several different lattice sizes (L = 
32fc, 64fc, 128fc, and 256fc) to perform a scaling analysis 
and to find the jamming concentration pj and the perco¬ 
lation threshold Pc at the thermodynamic limit (L —^ oo) 
(see, e.g. [2) 

PjA^) -PjA°°)) ^ 

where Vj and are the universal critical exponents for 
the jamming and percolation, respectively. In two di¬ 
mensions, = 1 [l| and iZc = 4/3 (2. 

The examples of scaling are shown in Fig. (jam¬ 
ming concentration, pj) and Fig.j^D (percolation concen¬ 
tration, Pc)- 


III. RESULTS AND DISCUSSION 
A. LKd model 

An example of a jammed state of fc-mers with defects 
on a square lattice with defects is shown in Fig.[TJi. Even 
at jamming filling, the percolation cluster looks rather 
sparse and lacy (Fig. [T]d). 

The results of the simulations are presented in Fig. |4] 
as phase diagrams on the plane {di,dk)- For any given 
fc, percolation can occur below the di{dk) curve. Even a 
small concentration of defects inhibits the deposition of 























(a) d = 7/9 ( 6 ) d= 11/18 


FIG. 2. Jamming states for different fractions of forbidden contacts (Cd model). The lattice size is 1024 x 1024, fc = 8, Zm = 18, 
fragments of the lattices with 128 x 128 sites are shown. Percolation is absent for large numbers of forbidden contacts, a > 14, 
(a) but is present at smaller numbers of forbidden contacts, z > 11, d = 11/18 (b). Online: Horizontal fc-mers are shown in 
red, vertical fc-mers are shown in blue, fc-mers of the percolation cluster are shown in gray, empty sites are shown in white. 
Print: gray-scale. 




FIG. 3. Examples of scaling dependencies for the jamming pj 
(a) and percolation (b) concentrations at different values 
of r = 0.5(1), 0.55(2), and 0.6(3). Here, = 1 and i/c = 4/3. 
k = 8. Cd model. The statistical error is smaller than the 
marker size when not explicitly shown. 



FIG. 4. Phase diagram presented as the concentration of the 
defects on the lattice, di, versus the concentration of the de¬ 
fects on the fc-mers, dk- The areas below the curves corre¬ 
spond to the percolation states. 


long objects. The curves separating percolating and non¬ 
percolating states are convex for short objects {k < 4) 
and concave for the larger objects (fc > 4). As the value 
of k increases, the area of the percolating region on the 
phase plane {di,dk) decreases faster than the product of 
0.5didk (Fig. [5]). We can expect that in the case of long 
objects, the area will be extremely narrow and located 
along the axes. The result suggests that for large objects, 
even a very small concentration of defects can inhibit the 
percolation when both kinds of defects are present. 

The area of the percolating state A in the phase plane 


(Figure [5]) may be well fitted by the function 

A = bQ + bik ^ + b2k + b^k ( 1 ) 

where bo = 0.1073 ± 0.0005, bi = 0.186 ± 0.002, 63 = 
0.111 ± 0.002, 63 = 0.0224 ± 0.0006, and = 0.9999. 
The limit of A = 0 corresponds to the critical length of 
kmax = 74 for the absence of percolation. This small 
value is unrealistic and is a result of an extrapolation 
based on data for rather short objects. Such extrapola¬ 
tion cannot be treated as reliable. 

The LKd model looks very promising to describe the 
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FIG. 5. Area of the percolation state on the phase diagram as 
a function of the inverse of the length of the fc-mers. The cir¬ 
cles correspond to the simulation data, the curve corresponds 
to a fitting 0 


real systems, nevertheless, it cannot give a sufficient es¬ 
timation of the critical length of the fc-mers at which 
percolation is impossible. 


B. Cd model 

Figure [5] presents the jamming concentration pj ver¬ 
sus the fraction of forbidden NN contacts, d, for different 
values of fc-mers. The jamming concentration of fc-mers 
decreases when the value of d increases. It is interesting 
that for fc > 2 and d = 1 (i.e. situation where NN is 
completely excluded) the value of pj is practically inde¬ 
pendent of k. At larger values of d the jamming concen¬ 
tration Pj increases with k. 

Figure [6] presents examples of the jamming pj (closed 
symbols) and percolation p,. (open symbols) concentra¬ 
tions versus the fractions of forbidden contacts, d, for 
k = 2 and fc = 16. Percolation was observed above a 
critical value of d, i.e., at d < dc- Increase of d up to dc 
resulted in an increase of the percolation concentration 
Pc- Percolation was confined between the lines Pj{r) and 
Pc{r) while at concentrations above pj percolation was 
suppressed by jamming. 

Figure [7] shows the percolation phase diagram pre¬ 
sented as the fractions of forbidden contacts, d, versus 
the length of the fc-mers. The critical value of dc goes 
through a maximum at fc « 4 as the value of k increases. 
The asymptotic behavior dc when fc —>■ oo may be fit¬ 
ted by a logarithmic function dc oc log k to the limiting 
length kmax ~ 12000 obtained by extrapolation of the ra¬ 
tio Pc/Pj 0 and kmax ~ 6000 obtained for models which 
consider the defects only on a lattice and only in the k- 
mers [s^. Thus the data presented in the present work 
confirms the previous estimations. Note that the loss of 



FIG. 6. Examples of the jamming pj (closed symbols) and 
percolation pc (open symbols) concentrations versus the frac¬ 
tions of forbidden contacts, d, for fc = 2 and k = 16. Percola¬ 
tion was observed above a critical value of d, i.e., at d < dc. 
The statistical error is smaller than the marker size. 



FIG. 7. Percolation phase diagram presented as the fractions 
of forbidden contacts d versus the length of the fc-mers. The 
filled squares correspond to the critical values of dc below 
which percolation is present. The dashed lines correspond 
to the logarithmic extrapolation of djk) dependence to the 
previous estimations kmax « 12000 Q, kmax ~ 6000 [30)| . 
The statistical error is smaller than the marker size. 


percolation for very long fc-mers is closely analogous to 
similar behavior observed for k x k squares [s^ . 

Figure m shows examples of percolation diagrams pre¬ 
sented as filling fraction p versus the length of fc-mer. For 
a fixed fraction of forbidden contacts, d, the percolation 
region is limited by the hatched area between the curves 
Pj{k) and Pc{k). The jamming concentration pj contin¬ 
uously decreases with increasing values of k whereas the 
percolation concentration pc goes through a minimum. 
The observed behavior is quite similar to that observed 
for d = 0 0113 . For the analyzed problem of cooper¬ 
ative sequential adsorption with a restricted number of 
lateral near-neighbor contacts, percolation was observed 
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in the interval kmin ^ k < kmax- The estimation gives 
kmin — kjYidx — S at d — 0.6, and k^i^ — 1, k^o^x — 
at d = 0.5. The inset to Fig. [8] shows the dependence of 
kmax versus d. At large values of kmax {kmax > 10), this 
dependence may be well approximated by the following 
logarithmic function 

^og{kmax) = a + bd, 

where a = 4.04±0.22, b = —4.93±0.57 and the coefficient 
of determination is = 0.9616. 


IV. CONCLUSION 

The jamming and percolation of linear fc-mers on a 
square lattice were investigated for two models. These 
models attempt to mimic the processes of deposition of 
elongated objects on real substrates. The first model KL^ 
deals with the random sequential adsorption of extended 
inhomogeneous objects onto a substrate with preliminary 
deposited impurities. These impurities inhibit the depo¬ 
sition of the objects. The inhomogeneity of the objects 
means that some sites are treated as insulating. Even for 
not very long objects, the presence of defects of two dif¬ 
ferent kinds (defects on the substrate and defects on the 
objects) at low concentrations prevents the percolation. 
The simulation suggests that for sufficiently large linear 
A:-mers percolation cannot occur even on an ideal lattice. 


The second model Cd is based on the cooperative se¬ 
quential adsorption model. This model takes into consid¬ 
eration interactions between the adsorbed objects. For 
this model the pj and Pc dependencies are controlled 
by the fractions of forbidden contacts d. The previ¬ 
ous works investigating the case of d = 0 (i.e., when 
all contacts are allowed) conjectured that the percola¬ 
tion of fc-mers is impossible if their length exceeds some 
critical value, kmax However, the estimations 

of the value of kmax used extrapolation to very long fc- 
mers (fc > 10^), as the direct simulation for such large 
objects was a very time-consuming task and was not 
very precise. The present work for the more general 
model at d > 0 estimates kmax more precisely for dif¬ 
ferent values of d. For the case of d = 0 the data from 
the present work confirms the previous estimations of 
kmax ss 12000 Q and kmax ^ 6000 [s^. The value of 
kmax decreases with increasing values of d. The log¬ 
arithmic dependence of the type log(fcmaa;) = a + bd 
{kmax > 0), where a = 4.03 ± 0.96, b = —4.93 ± 0.57, 
is obtained (i?^ = 0.9616). 
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